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Abstract. The temperature of the low-density intergalactic medium is set by 
the balance between adiabatic cooling resulting from the expansion of the universe, 
and photo-heating by the UV-background. A sudden injection of entropy from the 
reionization will increase the temperature of the gas, leading to a broadening of 
the hydrogen Lyman-a absorption lines produced in the IGM, and observed in the 
spectra of background quasars. We present a method based on wavelets to charac- 
terise objectively the line widths of such absorption lines. We use high resolution 
hydrodynamical simulations to demonstrate that the algorithm can detect changes in 
temperature of order of 50 per cent on scales > 5000 km s~^. We apply the method 
to a UVES/VLT spectrum of quasar 0055-269 (zcm = 3.7) and detect at the 99 per 
cent confidence level a sudden increase in temperature below redshift z ~ 3.3, which 
we interpret as evidence for Hell reionization. 



1 Introduction 

The intergalactic medium (IGM) can be observed in the spectra of distant 
quasars through resonant absorption in the hydrogen Lyman-a transition [Q. 
The absence of strong absorption throughout the spectrum (Gunn-Peterson 
trough) indicates that the IGM is very highly ionised, pm/pn ~ 10~^ below 
redshifts z ^ 6, but the nature of the first ionising sources remains uncertain. 

If the sources responsible for reionizing Hi are sufficiently hard, then Hen will 
be reionized at the same time. However, if for example galaxies were the 
main contributor to the UV-background, then Hen reionization may be sig- 
nificantly delayed until a population of harder sources appears. The large 
fluctuations in the Hen optical depth observed with STIS and FUSE Q sug- 
gest Hen reionization occurs late, around z ~ 3. There is supporting evidence 
from the hardening of the UV-background, as deduced from metal line ratios 
although these results remain controversial 

*Based on the data obtained from the OPC program 65.O-0296A (P.I. S. D'Odorico) at 
the VLT/Kueyen telescope, ESO, Paranal, Chile. 



An independent approach is to study the temperature of the IGM. Photo- 
heating and adiabatic expansion introduce a tight density-temperature relation 
[Q, T — Tq{p/ {p)y~^ , for the low density IGM responsible for producing the 
Lyman-a forest. A sudden injection of entropy resulting from reionization will 
increase Tq and make the gas nearly isothermal, 7 ~ 1. Detecting reionization 
through a sudden entropy increase has the advantage that one determines the 
epoch at which a significant fraction of the volume of the universe is overrun 
by the ionization front. 

Schaye et al. used hydrodynamical simulations to demonstrate that the 
density-temperature relation introduces a cut-off in the scatter plot of column 
density iVjjj versus line width b . They used a set of many high resolution 
simulations to calibrate the relation between the position of the cut-off and 
(To, 7). Applying the calibration to the cut-off measured in ten high-resolution 
spectra, they found evidence for a rise in To around a redshift z ~ 3.3, and 
an associated dip in 7, which they associated with Hen reionization. Ricotti, 
Gnedin & ShuU (|6|) applied a similar technique, but calibrated with pseudo 
hydrodynamical simulations, to published line lists, and found a similar tem- 
perature increase, although their error bars are large and their result is consis- 
tent with a non-evolving Tg as well. Bryan & Machacek also found evidence 
for a high value of Tq, but McDonald et al. Q did not find an increase in To 
around z — 3.3, although they used largely the same data as Schaye et al. 
Zaldarriage |^ applied a wavelet analysis similar to the one discussed here to 
look for temperature changes in the spectrum of QSO 1422+231 which might 
be a relic from reionization, and constrained them to be smaller than a factor 
of 2.5. 

2 Wavelet analysis 

A discrete wavelet is a localised function with a finite bandwidth. This makes 
wavelets useful for characterising line widths in a spectrum, since the amplitude 
of the wavelet will be related to the width of the line, and the position of the 
wavelet to the position of the line. The decomposition is moreover unique, for 
a given wavelet basis. 

Theuns & Zaroubi [|lOj used the Daubechies 20 wavelet to characterise 
temperature fluctuations in simulated Lyman-a spectra. They demonstrated 
how the wavelet amplitudes are large when the gas is cold and the lines narrow, 
and vice versa. They also showed how the cumulative distribution of wavelet 
amplitudes can be used to characterise the gas temperature, and to judge 
whether two stretches of spectrum have different temperatures or not. Here, 
we have used A, the maximum difference between the cumulative distribution 
of a stretch of spectrum, and the distribution for the spectrum as a whole, as 
a measure of how much Tq for a stretch, differs from {Tq). 

We have extended this method to attach a statistical significance to differ- 
ences in cumulative distributions. Given a line list, we construct mock spectra 
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Figure 1: Application of tlie wavelet method to a simulated spectrum (top 
panel) and QSO 0055-269 (bottom panel). In the simulations, the low redshift 
half is 50 per cent hotter than the high redshift half. This jump is detected at 
the 99.5 per cent level. A similar jump is seen in the observed spectrum. 



by scrambling the lines, thereby producing a spectrum where the lines have 
the same shapes as in the original spectrum, but any correlation in line widths 
that might result from temperature variations along the spectrum has been re- 
moved. We can now attach a statistical level of significance for the correlation 
between line widths - and hence the temperature of the gas - by comparing 
the statistics of the original spectrum with the scrambled one. 

In our simulation, we chose the photo-ionization rates for Hydrogen and He- 
lium such that Hi and Hei reionize around redshift z ^ 7 but Hen reionization 
is delayed until z ^ 3.4. In order to investigate the effect of different tem- 
peratures, we impose a power law equation of state, T ~ Tq {p/ {p))''^^ , 
on the simulation output before computing mock spectra. We have used 
(To, 7) = (1.5 X lO^K,^/?,) and (2.2 x IQ'^K,^/?,). These spectra are made to 
look like real data as much as possible, by adding photon and read-out noise, 
and broadening the lines with the instrumental profile of the spectrograph. 
Voigt profiles are fitted to the absorption lines with the same programme as 
used to analyse the data. 



3 Results 



Our results are shown in Fig. (|i|). The dot-dashed line (right hand scale) is the 
maximum difference between the cumulative distribution of wavelet coefficients 



for a window of size 8000 km at redshift z, and the spectrum as a whole. 
Histogram (left hand scale) is the significance (in %) for unusually large values 
of |A|, obtained from comparison with spectra constructed from randomised 
line lists. Negative values indicate gas colder than average, and vice versa for 
positive values. In the simulated spectrum (top panel), the first and second 
halves of the spectrum have Tq = 2.2 x IC^K and 1.5 x lO^K respectively. 
The wavelet statistic A detects this imposed jump, as can be seen from the 
sudden change in the sign of A. The significance of the change is computed 
to be > 99.5 per cent. The bottom panel applies the analysis to the spectrum 
of QSO 0055-269. We detect a sudden change in the value of A around 
z ~ 3.3, significant at the > 99.5 per cent. Note that no simulations are used 
to compute the level of significance: the method described here is solely based 
on the line list. 

We note that the high redshift gas is colder than the lower redshift gas. This 
is opposite to what one would expect from photo-heating in the optically thin 
limit, in which case the IGM will smoothly cool down. Therefore we suggest 
this sudden temperature increase is due to Hen reionization at redshift ^ 3.3. 
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